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Abstract 

We consider the three- loop singlet diagrams induced by axial- vector, scalar 
and pseudo-scalar currents. Expansions for small and large external mo- 
mentum q are presented. They are used in combination with conformal 
mapping and Pade approximations in order to arrive at results for the po- 
larization functions valid for all q 2 . Results are presented for the imaginary 
parts which are directly related to physical quantities like the production 
of top quarks or the decay of scalar or pseudo-scalar Higgs bosons. 

PACS numbers: 12.38.-t, 12.38.Bx, 14.65.Ha 



I. INTRODUCTION 

In the last few years a lot of effort has been devoted to developments of techniques 
which allow the evaluation of higher order corrections. Of special interest is thereby 
the evaluation of QCD corrections to two-point current correlators. Their knowledge 
immediately leads to a variety of important observables like the cross section <r(e + e~ — > 
hadrons) mediated by a photon or a Z boson or the decay of a scalar or pseudo-scalar 
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Higgs boson. Whereas at the one- and two-loop level exact results are known (for a 
review see IJ), until recently at 0(a 2 ) only expansions for large external momentum q 
respectively small quark mass m were available. 

In an approach was developed which leads to semi-analytical results for the three- 
loop polarization functions. The essence of this procedure amounts to the combination 
of the low- and high-energy analytical data for a polarization function through the use of 
the conformal mapping and Pade approximation suggested in |§-||. 

In a first step it was applied to the non-singlet diagrams induced by external vector, 
axial- vector, scalar and pseudo-scalar currents ||. In this paper this procedure will be 
applied to the corresponding singlet diagrams. They are often also referred to as double- 
triangle diagrams as the external currents are not connected through the same fermion 
line. This completes the knowledge of the three-loop current correlators at 0(a 2 ). Thus 
also the full mass dependence for the inclusive cross sections a(e + e~ — > tt) and the decays 
of a scalar or pseudo-scalar Higgs boson into quarks is available at this order. 

The method used in || heavily relies on the fact that the lowest particle threshold for 
the non-singlet graphs starts at q 2 = 4m 2 . In contrast to that, the singlet diagrams contain 
massless cuts. The solution of this problem is described in detail for the axial-vector 
correlator in Section ||. Also the essential ingredients are listed and the approximation 
procedure is briefly reviewed. Section |TTT| describes the treatment of the scalar and pseudo- 
scalar diagrams and finally results are presented in Section [TV]. Note that for the vector 
correlator there are no singlet diagrams at three-loop level according to Furry's theorem. 



II. SINGLET AXIAL- VECTOR CORRELATOR 



In this section the ingredients and the procedure for the construction of the Pade 
approximants of the axial- vector polarization function are presented in detail. The scalar 
and pseudo-scalar singlet diagrams are discussed in Section |TJ . 

Let us start with some definitions. It is convenient to introduce the dimensionless 
quantities 



Z = 4^' r = 4^' X = -TS' " = VY ^' (1) 

where a/s is the center of mass energy and m is the pole mass of the produced quark. 
x is a convenient variable in the high energy region and v represents the velocity of the 
quark. The axial-vector polarization function is defined through 

(-?V + <Wk) nV) +q»q„Tll(q 2 ) =i J dx e^(0\Tj;(x)f u (0)\0), (2) 

with j° = ip'Ju'j^ip. Only the transversal part IT a (g 2 ) will be considered in the following. 
In D = 4 — 2e space-time dimensions care has to be taken concerning the treatment of 75 
- especially in connection with the singlet diagrams. We follow the treatment introduced 
in and refer for more details to ||. 
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FIG. 1. Singlet or double-triangle diagrams. In the fermion lines either the quark ip or its 
isospin partner % may be present. 



In order for the axial anomaly to cancel, one has to take both members of a weak isospin 
doublet into account. It is therefore convenient to replace the current in Eq. (H) by 
is,/i = ^7^75^ — X7a*75X; where if) and x are isospin partners. The diagrams contributing 
to the singlet part, Il^(g 2 ), are depicted in Fig. [1] where in the fermion triangles either if> 
or x ma y De present. Note that for a degenerate quark doublet H|(g 2 ) vanishes. Having 
in mind the physical case (if>,x) = {t,b), however, we set = m and m x = in the 
subsequent analysis. H^g 2 ) is conveniently written in the form 

n^(g 2 ) = c F r(^) 2 4 2) -(g 2 ), (3) 

where Cp = (N 2 — 1)/ (2N C ) , and T is the trace normalization for an SU (N c ) gauge group. 
For QCD N c = 3 and T = 1 /2. 
The imaginary part of IT| (g 2 ), 

R a s (s) = 12vrlmn^(g 2 = s + ie), (4) 

normalized in analogy to the vector case, enters, e.g., the total inclusive cross section for 
the production of top quarks. However, also imaginary parts arising from the massless 
quark x contribute to Rs(s). The purely gluonic cut is zero according to the Landau- 
Yang- Theorem p. Let -Rsfc' a be the contribution of these massless cuts to Rg. Then the 
quantity 

rW) = j^-J 1 dr^lM (5 ) 
12tt z Jo r — z 

defines a function whose imaginary part for q 2 < Am? coincides with that of H^' a (q 2 ) 
and is zero for q 2 > 4m 2 . Its evaluation will be described below. With the help of this 
function we choose the overall renormalization condition for n^' a 



lim fn 



( 2 ).V)_n^V)) = . ( 6 ) 



In this difference the ln(— z) terms cancel for z — > which makes it possible to demand 
this QED-like renormalization condition. Alternatively one could use the MS scheme 
where only the 1/e poles are subtracted. 
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An important ingredient for the Pade approximation is the expansion of U^' a (q 2 ) 
for small external momentum. It is obtained by applying the so-called hard mass proce- 



dure [10 1 which provides a consistent expansion in q /m 



4 2) 'V) 



167T 



(2),a z n 



n>0 



(7) 



The first eight coefficients have been evaluated: 
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(2),tt 
5,0 



r »(2),a 



C 



c 



(2),a 
5,2 

(2),a 
5,3 



^(2),a 
°S,4 



c 



c 



(2),a 
S,5 

(2),a 



c 



(2),a 
5,7 



7>-a , ^/ _ r2 

471631 5215 56 

311040 + 13824 ^ 3 ~ 81 qm 
1178417 18179 _ 14 

8064000 + 92160 ~~ 75 qm 



15527720419 



1811719 13312 , 

+ „„ n &-^=L 



5120962560000 ' 14745600 ^ 165375 qm 
23054974995287 5965963 , 7696 , 

H C3 

178605 



811160469504000 70778880 



735623850793897673 



16260413 _ 

+ — c 3 



6600448 



in 



23817834479222784000 ' 264241152 ,u 252130725 qm 
1301167265336208772211 425212357 , 3007168 , 

+ Cs - I 



46002446022727434240000 9059696640 



173918745 



qm 



1 1 17124902949261 1725205473 



193004110999 



450455951454547036078080000 5218385264640 



(3 



132788224 
10956880935 



qm ; 



with £3 m 1.202056903 and l qm = \n(—q 2 /m 2 ). K a is a constant whose numerical value 
will be given below. 

A stringent constraint both for the real and imaginary part of Yi^ ,a {q 2 ) is set by the 
expansion for large external momentum. Recently, the large momentum procedure has 
been applied leading to the result M 1 : 



4 2) -v 



16tt 2 



K a + 



+ 



m 



2\ 2 



/ m 2 \ 3 
( m 2 \ 



185 385 m 2 
288 



16 



80 320 
yC 3 + — Cs_ 



380 



64C 3 + 



/296 

It ~ 



32C 3 U + 24/ 



3271 
243 



416 /280 \ 410 2 176 3 

"9" + \~2n + / m + ~27~ gm 27" 9m 



L In |8fl the results are listed in the MS scheme. 
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5584 



27 



105441373 



1340 

g ''qm 

2420 



c 3 



1660 



"qm 



101250 3 
6044237 _„ „ \ , 1177331 



112 Cs km + 



1350 



qm 



15542 3 
135 gm 



+ . . 



(9) 



with £3 defined above and ~ 1.036927755. 

The logarithms in Eq. (§) are due to the massless cuts and do not appear in the non- 
singlet diagrams 0. They spoil the procedure for constructing the Pade approximants 
developed in ||. Therefore, instead of dealing with the full polarization function, one 
may use Eq. (pj) and consider the following quantity: 



n 



(2),o 
S.mod 



(2),a/ 



(10) 



According to the definition of fl^' a (q 2 ) the ln(— z)-terms from the low energy expansion of 
Il^' a (g 2 ) are exactly canceled. Above z = 1, the imaginary part on the r.h.s. of Eq. (|l~0|) 



is determined by U s a (q 2 ) alone. This means that for z > 1, ^Smodil 2 ) contains all 
possible cuts of the double-triangle diagrams and one should subtract the massless ones 
by using Kg • 
quarks. 



s) in the region s > Am to get the production cross section for massive 

An analytic formula for Kgl' a (s) is available |ll|0- Nevertheless it is not possible 
to solve the dispersion integral in Eq. ( ]TD| ) analytically. On the other hand a purely 
numerical integration is excluded as the result contains ln(— z) terms for z — > which 
makes an expansion of the integrand with subsequent integration impossible. Let us 
therefore briefly describe the method we used for evaluation of Eq. (||). 
One may write Rfb^i 3 ) as ( reca U Eq. (|l|)) 



R 



(2),a, 
Sb \ 



3 - ln(4r) + R Sb 



(2),a, 



Then i?^' a (s) has a very simple limiting behaviour. For s — > it reads 

37 14 7 



1 r H r + . . . 

8 27 50 



'111 



(12) 



where the dots represent higher orders in r. The expansion of R Sb 
leads to 



(2),a, 



s) around = 2m 



In [11] a different kinematical region was considered. However, we continued the result to the 



region under consideration by using the translation table given in the appendix of [11] 
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4t» 



19 3 7 /l 7 

4 



■HVT^r) (27T--7rln2) +... 



(13) 



i.e., in this limit i?^' a (s) is a series in — r. Provided with this information we split 
the integral in Eq. $h into three parts: 



dr = /■ dr + r rfr + /■ dr m^i (14) 

r — z Jo r — z Jo r — z Js r — z 

and replace in a second step R^' a (s) in the interval [0, 5] by the expansion in fll2|) , in the 
interval [5, 1] by the one in flHf ). It turns out that the inclusion of the first 100 terms in 
the small energy expansion and the first 40 terms in the expansion around r = 1 leads 
to stable results in the range 5 = 0.65 . . . 0.80 with an accuracy of 13 to 14 digits. We 
will not quote numbers for the full flg' a (q 2 ) but only for the constant K a appearing in 
Eqs. © and (|): 

K a = -9.08040684374401 .... (15) 

The third kinematic region to be used for the Pade procedure is the threshold for the 
production of two massive quarks, z — > 1. In this region II^ (g 2 ) gets contributions from 
two sources: the cuts involving massive quarks and R^' a (s). It is strongly expected that 
the former starts at least with a term proportional to v in analogy to the non-singlet axial- 
vector correlator which follows the P-wave scattering solution of the Coulomb potential. 
Rgl' a (s) on the other hand has a smooth behaviour for s — ► 4m 2 . For z < 1, n^^g 2 ) 
is constructed in such a way that its imaginary part vanishes. However, the leading 
contribution of Yl ( g^ od (q 2 ) for z — > 1 + is given by 

r (2),a , 2 , 3 , / 1 \ / 19 . , „ _ 7 



n^ 2 ) = 7^ ln ( v T^J( v -y + 41n2 - 2 ^ + 6C3j + ... , (16) 

where the ellipses represent sub- leading terms in (1 — z). 

The construction of the Pade approximations divides naturally into four steps (for 
more detail we refer to 0): First, the threshold contribution has to be subtracted in all 
kinematical regions in order to have a polarization function that has a vanishing imaginary 
part for z — ► 1. Then, a new polarization function, n^£ od (q 2 ), is constructed whose high 
energy expansion contains no logarithmic terms any more. This must be done carefully 
in order not to destroy the behavior for z — > and z — > 1. In a third step the conformal 
mapping || 

(17) 



[l + u 



is used to transform the q 2 plane into the interior of the unit circle. Finally, a Pade 
improvement is performed in the new variable, u. 
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In only the constant and the m 2 /q 2 corrections in the high energy expansion 
have been included into the analysis. For the singlet diagrams meanwhile terms up to 
0((m 2 /q 2 ) 6 ) are available. This makes it necessary to modify the definition of the func- 
tion P{uj) for which the Pade approximation is performed. The natural extension of the 
definition given in 0] reads: 



{Au) n 1 / ~ (2)<a 2 ™ 1 1 / d? ~ (2) , a 2N 
n{ ' ~ (1+a;) 2 " I s > mod[q ' e J [WIW S ' mod[q ' 



+ u 



(4a;) 



where the index n > 1 indicates that the mass corrections of order (m 2 /q 2 ) n are included. 
The Pade approximants 

r.,. u \ a + aiu;H h a^o/ 

1 + Oiu; + • • • + fyu^ 

are then constructed from P n (— 1) and P^(0), (k = 0, 1, . . . , n + n — 1), where n is the 
number of moments (see Eq. (^)) used for the construction of the Pade approximation 
and P( k \0) = -rpfP n (ci>) | w =o- Taking into account all available information, i.e. n = 6 and 
n = 7, it is possible to construct approximants like [7/6], [6/7] or [8/5]. However, it turns 
out that the construction of P^(0) for large values of k suffers from huge cancellations. 
It is therefore necessary to evaluate the expressions in Eq. ([TJ|) with highest possible 
accuracy in order to arrive at reliable results for high-order Pade approximations. After 
all, the above mentioned 13 to 14 digits are enough to get stable results. Of course, 
also lower order Pade's have been evaluated both for consistency checks and to examine 
the convergence properties. We should mention that some of the Pade approximants 
develop poles for \u\ < 1 which result in poles in the physical z plane. Since this is not 
acceptable, only Pade approximants free from poles in the physical region are considered 
in the discussion of Section 0. We refrain from listing explicit formulae for U^' a (q 2 ) at 



this point and instead present in Section results for the imaginary part, R^' a ' : 



III. SCALAR AND PSEUDO-SCALAR CASE 

It is now straightforward to extend the procedure described above to the scalar and 
pseudo-scalar case. Here, in contrast to the singlet axial-vector contribution only the 
diagram with two massive triangles contributes. However, the cut through the two gluons 
does not vanish as it was the case for the axial-vector coupling, so that again there is a 
cut starting at z = 0. 

The polarization functions are defined through (k stands for s and p, denoting the 
scalar and pseudo-scalar case, respectively): 

q 2 U K (q 2 ) = i J dxe^(0\Tf(x)f(0)\0), (20) 
R K (s) = 87rimn K (g 2 = s + ie), (21) 
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where the currents are given by j s = ipip and j p = iipj^tp. As for j p one again has to deal 
with 75 in D ^ 4 dimensions, we again adopt the definition of 0, referring for details 
to[0. 

In analogy to (H) we define: 



n*(g 2 



In 2 jo 



1 , R { S' K (s) 

ar — — 



r — z 



[22] 



R(gJ' K (s) corresponds to the two gluon cut actually describing the Born decay of a scalar 
or pseudo-scalar Higgs boson to gluons [0,0: 



zr \ r 



gg 



(23) 



with 



f(r) 



arcsin 2 (y / r) 



log 



l+y'l-l/r 



i 2 



— 27T 



r < 1 
r > 1 . 



(24) 



Although these functions are quite simple an analytic integration is hard to perform. 
So we adopt the same procedure as for the axial- vector case and expand R^' K (s) for 
s — ► 0, where we take 200 terms into account, and for s — > 4m 2 , where 50 terms are 
enough to get a precision of 17 to 18 digits in the interval 5 G [0.65, 0.80]. 

In analogy to Eq. (0) the overall renormalization condition reads: 



Km(n?>'V)-ftt 







(25) 



In this scheme, the low-energy expansion of rig (g 2 ) looks as follows: 



C 



(2),« 
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(2),s 
S,2 

(2),s 
5,3 
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5,4 

(2),» 
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(2),s 
5,6 

(2),» 



b S,7 



4609 



721 , 
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2880 1152 ' 9 
2719121 10871 28 

5806080 + 36864 ^ 3 ~~ 135 qm 
519513881 1330021 1543 , 



3483648000 
2460910303 



7372800 
50939 , 



57480192000 409600 ' 14175 
8958934229477 2526649 , 
+ Ca 



14175 qm 
904 7 

221416 



2929190584320000 ' 27525120 ,J 5457375 ?m 

37498822356303853 5991294557 5844896 

"2999491158343680000 + 84557168640 ~~ 212837625 gm 
112110439141686419569 245566743541 , 20750416 , 

+ tttt^^ Cs - /. 



6118961963021107200000 4348654387200 



1064188125 



qm 
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(26) 



h qm 



11520 4608 
2468869 3969 7 

2419200 + 5120 ~ 15 ?m 
920009009 974281 _ 328 

1393459200 + 1638400 945 qm 
5898858645227 753259 3832 , 

+ TT^TTT C3 - — Z. 



12875563008000 
5523023003231 



1572864 
11702895 



Cs 



14175 ' qm 
3776 



16531587072000 ' 29360128 ,J 17325 qm 
63171728144529503 205215857 4266896 

249957596528640000 + 603979776 ^ 3 ~ 23648625 qm 
481264894165689829721 513043585411 32459264 , 

Cs- z 
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qm j 



(27) 



where !!<? (g 2 ) and Cg„ are defined in analogy to Eqs. @ and (|7|), and 



K s 



0.62280338337755 



K p 



1.81359971877046 



(2f 



The high energy expansion terms up to 0(1/ z 4 ) are already listed in [12 in the 
MS scheme. We have added the (m 2 /s) 5 and (m 2 /s) 6 mass correction terms. In the 
renormalization scheme defined in Eq. fl25|) the result reads (l qfl = ln(— g 2 //i 2 ), with \x 
being the renormalization scale): 
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(30) 



It should be noted that in the scalar case the singlet polarization function explicitly 
depends on the renormalization scale fi. For the approximation procedure the choice 
fi 2 = m 2 will be adopted. The imaginary part, however, is independent of \i. 



The analogue of Eq. ([U]) reads 



ll S,mod\H 



n 



4 2) 'V 



(31) 



The leading threshold term of this function both for k = s and p will be determined 
in complete analogy to the axial- vector case (see Eq. (|16D ), despite the fact that for the 
non-singlet part in the pseudo-scalar case it originates from the S"-wave solution of the 
Coulomb potential. If by this procedure an essential contribution to the threshold part is 
missing, different Pade results should develop a large spread close to v = 0. One obtains 
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16vr s 
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(32) 



The approximation procedure is applied in complete analogy to the axial-vector cor- 



relator. As for the scalar and pseudo-scalar case also the eighth moment, Cg 
available, Pade's like [7/7], [8/6] or [6/8] may be evaluated. 
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FIG. 2. The imaginary parts Rg' a , Rjp' s and Rg' p of (a) the axial-vector, (b) scalar 
and (c) pseudo-scalar singlet diagrams, respectively. For the abscissa the variable x = 2m/ ^/s 
is chosen. Solid line: Pade result; wide dots, dashes and narrow dots: (m 2 /q 2 ) n -expansion for 
n = 0, n = 1, . . . , 5 and n = 6, respectively; dash-dotted line: purely massless cuts R^' a , 
R^gg' s and R g 2 g' p - (d) shows the difference between the solid and the dashed dotted line (i.e., 
the contribution of the massive quarks) of (a), (b) and (c) as solid, dashed and dotted line, 
respectively. 



IV. RESULTS 

In Fig. (a)-(c) the results for the imaginary part of n^' ) ' a , IIg and Hg p (solid 
lines), together with the first seven terms of the high energy expansion (dashed and 
dotted lines) are shown as functions of x = 2m/y/s. Recall that in the displayed region, 

(2) (2) 

< x < 1, ImlTg = \mI\. Smod . Therefore, if one is interested, e.g., in (inclusive) 
production of the heavy quarks only, the corresponding massless cuts (depicted as dash- 
dotted lines) have to be subtracted. The resulting curves are shown in Fig. (d). 

As was already pointed out in ||, the imaginary part of the axial- vector singlet con- 
tribution starts at 0(m 6 /s 3 ). Therefore, Rg' a is rather small below x = 1/2. On the 
other hand, above this value the four-particle threshold at y/s = Am is expected to hi- 



ll 



hibit convergence of the high energy expansion. However, there is still an agreement 
of the (m 2 /s) 6 -terms with the semi-analytical Pade result up to x ~ 0.7. Also in the 
scalar and pseudo-scalar case convergence of the high energy expansion is quite fast up 
to x ~ 0.7. . . 0.8. For the pseudo-scalar case even very close to threshold the difference 
to the Pade result is rather small. 

Minor differences among the various Pade approximations are visible only close to 
threshold. In Fig. || (a)-(c) for each correlator 30 Pade approximants with different input 
from q 2 — > and q 2 — > — oo are plotted as functions of the quark velocity v = \Jl — Am 2 / s. 
The dotted curve shows the (m 2 /s) 6 -terms of the high energy expansion for comparison. 
The spread of the different Pade approximations serves as a measure for the uncertainty 
of the procedure. The dash-dotted line again corresponds to the purely massless cuts 
which largely dominate the singlet contribution in all three cases. For the axial-vector 
and the scalar case the differences Rg' a — R^' a , respectively, R^' s — Rf) ,s are remarkably 
small below v — 0.5 which confirms our previous assumptions on the threshold behavior 
in these cases (see Fig. || (d)). As was expected the analogue quantity in the pseudo- 
scalar case grows steeper for v — » 0, but still the stability of the Pade result justifies the 
threshold assumption. We even modified artificially Eq. (|32| ) by ±10% and observed a 
drastic destabilisation of the Pade approximants. This also gives quite some confidence 
to the ansatz for U^' p at threshold. 

In the following handy approximation formulae for the imaginary parts R^'^R^ 
and R { s ] ' p are given. To keep the formulae as short as possible a Pade approximant 
containing the high energy input up to 0{m 2 /s) is selected. After subtraction of the 
known analytic results the (small) remainder, which vanishes both for v — and v — 1 is 
approximated using Legendre polynomials. The results read: 




+ 50p 3/2 (l -p) [-0.0789 PoO) - 0.0869 Pi (p) + 0.1099 P 2 (p)] , (33) 




+ 50p 3/2 (l - p) [-0.193 P„(p) + 0.111 Pi(p) + 0.0277 P 2 (p)] , (34) 
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FIG. 3. B§ ),a , R S ),S > R S ] ' P as functions of v = yjl-4m' 2 /s. Solid and dashed lines: 
30 different Pade curves; dotted: high energy expansion up to (m 2 /s) 6 ; dash-dotted: purely 
massless cuts, (d) is the analogue of Fig. || (d), plotted over v. 



+ 50p 3/2 (l - p) [1.586 P Q 0) - 1.822 P^p) + 0.525 P 2 (p)] , (35) 

with p = (1 — v)/(l + v). Note again that the corresponding massless cuts (see Eqs. (|TT|) 
and Q2"3|)) have to be subtracted in order to find the rates for, e.g., tt production. 

To summarize, the polarization functions for axial-vector, scalar and pseudo-scalar 
singlet current correlators have been computed. To this end the knowledge of n(g 2 ) for 
q 2 — > 0, q 2 — ► — oo and q 2 — > 4m 2 has been exploited in combination with conformal 
mapping and Pade approximation. Results have been presented for the imaginary parts 
and handy approximation formulae have been provided. This completes the evaluation 
of three-loop polarization functions at 0(a 2 ) as in || the non-singlet contributions were 
considered. 
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